> 

o 






X 



IFUP-TH 8/97 

On the evaluation of universal non-perturbative constants in 0{N) 

a models. 



Massimo Campostrini, Andrea Pelissetto, Paolo Rossi, and Ettore Vicari 
Dipartimento di Fisica dell'Universita and I. N.F.N. ,1-56126 Pisa, Italy 

(February 1, 2008) 



Abstract 



D 

We investigate the relation between on-shell and zero-momentum non- 
perturbative quantities entering the parametrization of the two-point Green's 
function of two-dimensional non-linear 0{N) a models. We present accurate 
estimates of ratios of mass-scales and renormalization constants, obtained by 
an analysis of the strong-coupling expansion of the two-point Green's function. 
^^ [ These ratios allow to connect the exact on-shell results of Refs. ]|J2[ with 

Q I typical zero-momentum lattice evaluations. Our results are supported by the 

t^^ ' 1/A^-expansion. 
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Introduction. Physical quantities which are independent of coordinates and carry no 
physical dimensions, like mass or amplitude ratios, are the best candidates for scheme- 
independent and/or numerical determinations in quantum and statistical field theories. 

In two-dimensional non-linear 0{N) a models some exact results concerning the on- 
shell (large-distance in the Euclidean space) behavior of the two-point spin-spin Green's 
function are known. Exact formulas have been presented for the on-shell mass-A-parameter 
ratio [Q], and for the constant Ai of the 0(3) a model 0, which is defined starting from a 
parametrization of the large-momentum asymptotic behaviour of the on-shell renormalized 
two-point correlation function. No exact off-shell results are known. In this letter we study 
the relation between on-shell and zero-momentum quantities related to the two-point Green's 
function. This will allow, by using the above-mentioned exact on-shell results, to determine 
very accurately (the error being of order 10~^) the asymptotic behavior (i.e. for f3 — >■ 
cx)) of the non-perturbative quantities which parametrize the behavior of the two-point 
Green's function at small momentum, such as the second-moment mass and the magnetic 
susceptibility. 

For the sake of generality, let us discuss the general 0{N) model on a square lattice with 
nearest-neighbor action 

S = -N(3Y.s{x)-s{x + ij), (1) 

X,fl 

where s{x) ■ s{x) = 1. We have introduced a rescaled inverse-temperature /?, and we shall 
use the short-hand 

^ - 2 
" = 27NP- (2) 

We consider the Fourier transform of the bare spin-spin two-point function 

G{p;(3) = Y.e'^^{s{x)-m)- (3) 

X 

In the continuum limit and in the large (euclidean) momentum regime a standard one-loop 
calculation gives 

Af-l 

here a{p) denotes the (one-loop) running coupling constant 

«(P) = T— ^] , (5) 

1 -I- a m ap 

where a is the lattice spacing. The corresponding renormalized Green function G{i{p,M) 
defined by an on-shell renormalization condition is obtained by requiring 



where M is the physical mass-gap. The renormahzed correlation function can then be 
parametrized in the large-momentum regime as 

GRiv^m^^Hp)-^^- (7) 

Notice that Eq. (|^) is the definition of Ai(A^). This constant depends crucially on the 
renormalization condition one adopts. In turn the bare two-point function (^ depends 
explicitly on the coupling. In the low-momentum (large-distance) regime and in the scaling 
region it can be parametrized by 

Renormalization-group arguments lead to the following expression for the (universal) asymp- 
totic behavior of the mass-gap: 



M{/3) = -R{N)a~N~2e-a[l + 0{a)], (9) 

and of the on-shell renormalization constant 

N-l 

Z{p) = C{N)aN-2[i + 0{a)]. (10) 

The constant R{N) is not universal and it can be easily computed from the exact result 
of Ref. m by calculating the appropriate A-parameter ratio, which can be obtained by a 
simple one-loop calculation. For the standard nearest-neighbor action on the square lattice 
one finds 

1 



R{N)=(-] " r(l + ^)" e2(^-2)v^. (11) 



The constant C{N) is universal, that is independent of the lattice regularization. By using 
Eqs. d), (0), (ID and d), it can be put in relation with Ai(A^): 

For A^ = 3, using the exact result of Ref. 0, i.e. 

^i(3) = ^, (13) 

we get 

C(3) = SttI (14) 

The main purpose of the present letter is to compute the relation between R{N) and 
C{N) and the corresponding zero-momentum quantities, which are much easier to compute 



in Monte Carlo lattice simulations. We will obtain results which are quite accurate, although 
not exact. 

Typical lattice calculations lead to estimates of moments of the two-point function: 

m2^^Y.i^'y{m-s{x)). (15) 

X 

In particular x = ^^o- If we now parametrize the function G{p; (3) around p = by 

{Mg is by definition the inverse of the so-called second- moment correlation length), we 
obtain the relationships 

Za = xMl = ^. (18) 

1712 

Standard renormalization-group arguments lead again to 

MG{(3) = RG{N)a~^^e'a[l + 0{a)], (19) 

N-l 

Zg^ = CG{N)aN-2 [1 + 0(a)] . (20) 

Cg{N) and RgIN) differ from the corresponding quantities C{N) and -R(A^) (except at 
N = oo where the theory is Gaussian). In order to calculate for these two quantities, we 
will thus investigate the dimensionless ratios 



a^O Mq oi-*0 Zg 



These quantities are not exactly known and we will provide here rather accurate strong- 
coupling estimates, supported by a 1/A^ analysis. 

Strong-coupling estimates of dimensionless RG invariant quantities. As shown 
in Ref. |^, strong-coupling analysis may provide quite accurate continuum- limit estimates 
when applied directly to dimensionless renormalization-group invariant ratios of physical 
quantities. The basic idea is that any dimensionless renormalization-group invariant quantity 
R{P) behaves, for sufficiently large /3, as 

R{(3) -R* ^ M{(3)\ (22) 

where R* is its fixed point (continuum) value and M{(3) goes to zero for /3 ^ oo. Hence 
a reasonable estimate of R* may be obtained at the values of /? corresponding to large but 
finite correlation lengths, where the function R{(3) flattens. This is essentially the same idea 
underlying Monte Carlo studies of asymptotically free theories, based on the identification of 
the so-called scaling region. Strong-coupling estimates of physical quantities may be obtained 



by evaluating approximants of their strong-coupling series at values of j3 corresponding to 
reasonably large correlation lengths, e.g. ^ ^ 10. Scaling is then checked by observing the 
stability of the results varying [3. 

In a strong-coupling analysis it is crucial to search for improved estimators of the quanti- 
ties at hand, because better estimators can greatly improve the stability of the extrapolation 
to the critical point. Our search for optimal estimators was guided by the large- A^ limit of 
lattice 0{N) a models which is a Gaussian theory. We chose estimators which are perfect 
for A^ = oo, i.e. do not present off-critical corrections to their critical value. 

On the lattice, in the absence of a strict rotation invariance, one may define different 
estimators of the mass-gap M having the same critical limit. On the square lattice one 
may consider ^ obtained by the long-distance behavior of the side wall-wall correlation 
constructed with G{x), or equivalently the solution of the equation G~^{i^, 0; /?) = 0. At a 
finite order q of the strong-coupling expansion, the wall-wall spin-spin correlation function 
Gw{z) at a distance larger than g/3 exponentiates exactly, i.e. for \z\ > q/3 it can be written 
as 

G^(^;/3) = A(/3)e-''('5)l^l. (23) 

In the context of a strong-coupling analysis, it is convenient to use another estimator of the 
mass-gap derived from /u(/?) 0: 

M^{(3) = 2 (cosh/i(/5) - 1) . (24) 

Moreover, by comparison with Gaussian model, we consider the following estimator of Z 

Zsil3) = 2A{/3)smh^i{(3). (25) 

In practice, when the strong-coupling expansion of G{x; 13) is known to order q, M'^{j3) can 
be determined up to about 2g/3 orders 0], and the same precision is therefore achieved in the 
determination of Zs{j3). Estimators of the zero- momentum mass Mq and renormalization 
constant Zq can be easily extracted from Eqs. (|1^ and (|^). So in order to estimate Sm 
and Sz one should study the continuum limit of the ratios 



Sm{(3) ^ ^, Sz{(3) ^ ^. (26) 

In the large- A^ limit Sm{P) = Sz{P) = 1 independently of /?. Quantities having the same 
properties of M^ and Zg can be conceived also on the honeycomb and triangular lattices 0, 
thus leading to analogous definitions of Sm{P) and Sz{(3). 

Analysis of the strong-coupling series. We have analyzed the strong-coupling series 
of Sm and Sz on the square lattice, where the available series |^ are of the form 1 -|- 
(5^ Ya=o (^iP^ in both cases. An analogous analysis has been performed on the honeycomb 
and triangular lattices within their nearest-neighbor formulations using the available series 
of the two-point function |0|. The analysis of strong-coupling series calculated on different 
lattices offers a possibility of testing universality. On the other side, once universality is 
assumed, it represents a further check for possible systematic errors, whose estimate is 



usually a difficult task in strong-coupling extrapolation methods such as those based on 
Pade approximants and their generalizations. 

In Table | we present some details of the results obtained on the square lattice. There 
we report estimates of Sm and Sz at various values of /3, where the correlation length is 
reasonably large. Such estimates of Sm and Sz have been obtained by resumming the strong- 
coupling series of {Sm — l) / (3^ and iSz — l) //5^ by Dlog-Pade approximants (DPA's), in 
which the standard Pade resummation is applied to the series of the logarithmic derivative, 
and then the original quantity is reconstructed. This method of resummation turned out to 
give the most stable results (we also tried simple Pade approximants and first order integral 
approximants). For a nth order series, we considered \l/m\ DPA's having 

/ + m + l>n-2, /,m>f-2. (27) 

As estimate at a given j3 we took the average of the values of the non-defective approximants 
constructed using all available terms of the series. As an indicative error we considered 
the square root of the variance around the estimate of the results from all non-defective 
approximants specified above. This quantity should give an idea of the spread of the results 
from different approximants. Approximants are considered defective when they present 
spurious singularities close to the real axis for Re/? ^ /3. 

The precision of the results is satisfactory even for values of (3 where the correlation 
length is quite large. Furthermore scaling is well verified. Then, assuming scaling, we 
extracted estimates of the corresponding continuum limit, which are reported in Table ||. 
There we report results obtained on the square, honeycomb and triangular lattices, and for 
several values of N [N = 3, 8, 16, where the last two large values of N has been considered 
in order to make a comparison with the large- iV analysis, see later). Errors represent a 
rough estimate of the uncertainty, which is quite small. Universality among different lattice 
formulations is well verified. Our final estimates for A^ = 3 are 

Sm = 0.9987(2), Sz = 1.0025(4). (28) 

There are some estimates of Sm obtained by high-statistics Monte Carlo simulations 
that are worth being mentioned for comparison. Monte Carlo simulations at A^ = 3 |^ 
gave Sm = 0.9988(16) at /3 = 1.7/3 = 0.5666... (^ ~ 35), and Sm = 0.9982(18) at (3 = 0.6 
(^ ~ 65), leading to the estimate Sm = 0.9985(12). From the data of Ref. [^ one derives 
Sm = 0.996(2) for A^ = 3 and Sm = 0.9978(8) for A^ = 8. These numbers compare well with 
our strong-coupling calculations, which appear to be much more precise. 

The two-point Green's function at small momentum. The fact that both Sm ~ 1 
and Sz ^ I should not come as a surprise. Indeed it was shown in Ref. that in the 
region p"^ ^ Mq the spin-spin two-point function is essentially Gaussian with very small 
corrections. The inverse two-point function can be studied around p^ = by expanding it 
in powers of p^: 



G-\p) = Z^'M, 



G 



,^P"^ f P" 
1 + —^ + > , C 



M'g t. ' \Ml, 



(29) 



Analysis based on various approaches have shown that in two- and three-dimensional 0{N) 
models the following relations hold [0,0,|[ 



C2 < 1, Q -C C2 for i > 2. (30) 

Then neglecting all Cj, z > 3 and terms of order c^ one may write the following expression 
for the inverse two-point function 

n2 



G-\p)^Za'M^ 



p2 



1 , P 



(31) 



which should give a good approximation of the two-point Green's function in the region 
Ip^I ^ Mq. As a consequence one obtains the following approximate relations 

Sm^I + C2, (32) 

Sz^l-2c2. (33) 

In order to have a check of Eqs. (|3^) and (^3|), we have estimated C2 by a strong-coupling 
analysis. Again guided by the large- iV limit, we considered the following estimator of C2 

In the continuum limit C2(/3) —>■ C2- In the large- A^ limit C2(/3) = C2 = independently of f3 
for all square, honeycomb and triangular lattices. We mention that on the square lattice, 
where G{x) has been calculated up to 21st order 0, the available series of C2 is of the form 
P'^ J2l=o (^iP^ ■ Details of the analysis of C2 on the square lattice can be found in Table |. 
Final estimates are reported in Table 0, where also results obtained on the honeycomb and 
triangular lattices are reported. Numerically in the case of 0(3) we obtained C2 — —0.0012, 
which nicely confirms Eqs. (|3^) and (p3|). 



1/A^ expansion. We have also evaluated Sm and Sz in the context of the 1/A^ expansion, 
which was found to be a fairly accurate approach to the evaluation of amplitude ratios in 
two-dimensional 0{N) models for N > 2. An analytic computation of the two-point function 
in the region p^ + M^ ^ leads to the following result: 






l + ^(log— -7£;-3)+0(l/Ar2^ 



(35) 



27ra 
and therefore 

XiiN) = 1 + ^ (log ^ + 7ir - 2) + 0{1/N'). (36) 

By evaluating the 0(1/A^) contribution to the on- shell- renormalized self-energy in the region 
around p = 0, one finds 

5m-1-^^T^^0(^). (37) 

...,,££HiZ21£,o(-l), (38) 

0.00619816 ^ / 1 \ ,,^, 

0.00023845 ^/ 1 \ 
-3 = ^^ + 0(^), (40) 

0.00001344 ^/ 1 \ ,,,, 



etc... The coefficients of the 0{1/N) terms show consistency with Eqs. (|32|) and (|33D within 
errors of order 10~^. Inclusion of C3 ~ 2.4 ■ 10~^/iV would squeeze the error to O(10~^). 
Strong-coupling estimates reported in Table y clearly approach the large-A^ asymptotic 
regime predicted by the above equations. In particular quantitative agreement (within the 
uncertainty of our strong-coupling calculations) is found at A^ = 16. This represents a further 
check of the analysis employed in order to get strong-coupling estimates in the continuum 
limit. 

It is worth mentioning that C2, Sm and Sz have been also calculated to 0{e^) within 
the (f)^ formulation of 0{N) models in 4 — e dimensions |^. The approximate relations ( |52D 
and (p3D are con&med even in the e-expansion, whose validity should not be related to the 
specific value of N. Furthermore a semi-quantitative comparison, inserting the value e = 2 
in the e-expansion formulae, provides the correct order of magnititude. 

Conclusions. Putting together the exact formulas ([TT|) and ([l^) and our strong-coupling 
estimates of the ratios Sm and Sz, we arrive at the following results 

Rg{3) = 4P = 80.139(8), (42) 

Cg{3) = C(3) xSz = 93.25(3). (43) 

For comparison we mention the existing Monte Carlo results concerning the ratio 
CG(3)/i?G(3)^ Ref. quotes CG{S)/RG{3f = 0.0146(11), which has been obtained by 
employing finite-size-scaling based techniques allowing to reach correlation lengths up to 
0(10^). Ref. § quotes CG(3)/i?G(3)2 = 0.0138(2), which has been obtained by standard 
Monte Carlo simulations up to ^ ~ 130, and where the error is just statistical. An attempt to 
estimate the systematic error due to violations of asymptotic scaling would give the number 
0.0138(2) (7)[|, where the second number within brackets is the systematic error estimated 
by us. These numbers compare very well with our corresponding estimate derived from 
Eqs. dg) and (H), i.e. CG(3)/i?G(3)^ = 0.01452(5). 



^ In Ref. B the estimate of the ratio CG(3)/i?G(3)^ has been obtained by using the Symanzik tree- 
improved action and considering the so-called energy scheme. In order to estimate the systematic 
error due to violations of asymptotic scaling, we have considered the difference between the results 
obtained by using two-loop and three-loop formulas in the fit of Monta Carlo data. Two-loop and 
three-loop formulas lead to 0.0145(2) and 0.0138(2) respectively. 
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TABLES 



TABLE L Estimates of Sm-, Sz and C2 obtained from DPA's of the strong-coupling series on 
the square lattice evaluated at various values of (3 where the correlation length ^ is reasonably large. 
For example at iV = 3: i{(3 = 0.45) ~ 8, ^{(3 = 0.50) ~ 11, ^{(3 = 0.55) ~ 25, and ^(/3 = 0.60) ~ 65; 
at iV = 8: C(/3 = 0.50) ~ 5, ^{(3 = 0.55) ~ 8, and i{f3 = 0.60) ~ 12. 



N 




[3 = 0.45 


/3 = 0.50 


(3 = 0.55 


13 = 0.60 


3 


{Sm - 1) X 103 
{Sz - 1) X 103 

C2 X 103 


-1.01(3) 

2.04(6) 

-1.16(4) 


-1.2(1) 

2.4(2) 
-1.3(2) 


-1.1(2) 

2.2(4) 

-1.4(3) 


-0.9(3) 

1.4(7) 
-1.3(4) 
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(Sm - 1) X 103 
{Sz - 1) X 103 

C2 X 103 


-0.54(1) 

1.11(3) 

-0.58(2) 


-0.59(2) 

1.2(1) 
-0.63(5) 


-0.57(5) 

1.2(1) 
-0.6(1) 


-0.5(1) 

1.1(2) 

-0.6(2) 


16 


(Sm - 1) X 103 
{Sz - 1) X 103 

C2 X 103 


-0.32(4) 

0.7(1) 
-0.31(1) 


-0.35(6) 

0.7(1) 
-0.34(3) 


-0.4(1) 
0.8(2) 
-0.35(6) 


-0.4(2) 
0.9(4) 
-0.35(10) 



TABLE IL We report estimates of Sm^ Sz and C2 from the strong-coupling expansion on the 
square, honeycomb and triangular lattice, and 1/A^ expansion of the continuum formulation of the 
non- linear 0{N) a model. Final strong-coupling estimates are taken at /3- values corresponding to 
^~ 10. 



N 




Sm 


Sz 


C2 


3 


square 


0.9988(2) 


1.0024(4) 


-1.3(2)xl0-3 




honeycomb 


0.9986(3) 


1.0027(4) 


-1.2(2)xl0-3 




triangular 


0.9985(5) 


1.003(1) 


-1.2(3) X 10-3 




0{1/N) 


0.9978 


1.0044 


-2.07x10-3 
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square 


0.99943(5) 


1.0012(1) 


-0.6(1) X 10-3 




honeycomb 


0.9994(1) 


1.0011(2) 


-0.7(1) X 10-3 




0{1/N) 


0.99919 


1.00164 


-0.77x10^3 


16 


square 


0.9996(1) 


1.0008(2) 


-0.35(5)xl0-3 




honeycomb 


0.9997(1) 


1.0006(1) 


-0.36(3) X 10-3 




0{1/N) 


0.99960 


1.00082 


-0.39x10-3 
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